Exerccios 4 - Anlises Exploratrias

4.2 Cervejas
Uma amostra de 30 estudantes foi indagada sobre seu tipo de cerveja preferida, com o seguinte resultado.
1. Represente este resultado como um grfico de barras e um dotplot (funo dotchart).
grfico de barras :
cervejas <-c("chope","lata","garrafa","chope","garrafa", "garrafa","lata","lata","nenhuma","lata","garrafa","garrafa", "garrafa","lata","lata","lata","garrafa","lata","chope","nenhuma", "garrafa","garrafa","garrafa","chope","garrafa","garrafa","chope","garrafa","lata","lata")
class(cervejas)
cervejas1<- factor(cervejas)
cervejas2<-table(cervejas1)
barplot(cervejas2)
barplot( cervejas2 , xlab="Tipo de cerveja preferida", horiz=F, las=1)


dotplot  :
# Grfico de barras por dotchart
par( omd=c(0,1,0,1) )
cervejas <-c("chope","lata","garrafa","chope","garrafa", "garrafa","lata","lata","nenhuma","lata","garrafa","garrafa", "garrafa","lata","lata","lata","garrafa","lata","chope","nenhuma", "garrafa","garrafa","garrafa","chope","garrafa","garrafa","chope","garrafa","lata","lata")
class(cervejas)
cervejas1<- factor(cervejas)
cervejas2<-table(cervejas1)
class(cervejas2)
cervejas3<- matrix(cervejas2)
colnames(cervejas3)<- c("Tipo de cerveja preferida")
rownames(cervejas3)<- c("chope", "garrafa","lata", "nenhuma")
dotchart(cervejas3)

2. Qual tem maior razo dado/tinta?
As duas figuras apresentam bem os dados.

4.3 Caixetais
1.  Construa um histograma do dap1) dos fustes dos caixetais.
caixeta<- read.table("caixeta.csv", header=T, sep= ",", dec=".")
head(caixeta)
caixeta2<-caixeta
caixeta2$dap<- caixeta$cap/pi
head(caixeta2)
hist( caixeta2$dap )



2.   Construa histogramas da altura das rvores para os diferentes caixetais ('local').
# Meu objeto  hist( caixeta2$h )
hist( caixeta2$h[ caixeta2$local == "chauas" ], xlab="Altura da rvore", ylab="Freqncia", main="Altura da rvore - Chaus", col = "blue")


hist( caixeta2$h[ caixeta2$local == "jureia" ], xlab="Altura da rvore", ylab="Freqncia", main="Altura da rvore - Jureia", col = "green")





hist( caixeta2$h[ caixeta2$local == "retiro" ], xlab="Altura da rvore", ylab="Freqncia", main="Altura da rvore - Retiro", col = "yellow")

3.   H diferenas entre as estruturas (distribuio de tamanhos) dos caixetais?
Sim, h diferenas sendo que a regio da Jureia se aproxima de uma distribuio normal com mdia em 100.





4.4 Eucaliptos

1.	Utilize o grfico boxplot para analisar o DAP de rvores de E. grandis em funo das variveis regies (regiao) e rotao (rotacao).


egrandis<- read.table("egrandis.csv", header=T, sep= ";", dec=".")

head(egrandis)

# Boxplot do Dap

boxplot(egrandis$dap)

# Boxplot do Dap em funo da regio e rotao

boxplot( dap ~ regiao*rotacao, data=egrandis)


















2.	Avalie a normalidade da altura do conjunto total de rvores com um grfico quantil-quantil contra a distribuio normal.


Com exceo dos valores das extremidades,  a varivel altura parece seguir uma distribuio normal.

4.5 Mais Caixetais
1.	Analise a relao dap-altura ('dap' e 'h') em funo do caixetal (local) com a funo plot, mas somente para as rvores 2) de caixeta (Tabebuia cassinoides).

coplot( h ~ dap | local, data=caixeta2)

coplot( h ~ dap | local, data=caixeta2 , panel= panel.smooth)





2.	Para a mesma relao do item anterior, verifique linearidade com a funo scatter.smooth

scatter.smooth( caixeta2$dap, caixeta2$h , col="red" )













3.	Utilizando o pacote lattice, analise a relao dap-altura ('dap' e 'h') em funo do caixetal (local), mas somente para as rvores 3) de caixeta (Tabebuia cassinoides).
# Pacote Lattice

library(lattice)
caixeta<- read.table("caixeta.csv", header=T, sep= ",", dec=".")
head(caixeta)
caixeta2<-caixeta
caixeta2$dap<- caixeta$cap/pi
head(caixeta2)
class(caixeta2)
xyplot( h ~ dap | local , data=caixeta2 )


